Local Charge Excesses in Metallic Alloys: a Local Field Coherent Potential 

Approximation Theory 
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Electronic structure calculations performed on very large supercells have shown that the local charge 
excesses in metallic alloys are related through simple linear relations to the local electrostatic field 
resulting from distribution of charges in the whole crystal. 

By including local external fields in the isomorphous Coherent Potential Approximation theory, 
we develop a novel theoretical scheme in which the local charge excesses for random alloys can 
be obtained as the responses to local external fields. Our model maintains all the computational 
advantages of an isomorphous theory but allows for full charge relaxation at the impurity sites. 
Through applications to CuPd and CuZn alloys, we find that, as a general rule, non linear charge 
rearrangements occur at the impurity site as a consequence of the complex phenomena related with 
the electronic screening of the external potential. This nothwithstanding, we observe that linear 
relations hold between charge excesses and external potentials, in quantitative agreement with the 
mentioned supercell calculations, and well beyond the limits of linearity for any other site property. 
PACS: 71.23.-k, 71.20.-b 
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I. INTRODUCTION 

For many years the Coherent Potential Approxima- 
tion (CPA) [2] has been widely used for calculating the 
electronic properties of random metallic alloys. There 
are many reasons for such a fortune: as it has been al- 
ready pointed out, the CPA allows for a very careful 
determination of spectral properties [3] and Fermi sur- 
faces [4]. However, for the arguments discussed in the 
present paper, two other aspects of the CPA are particu- 
larly amenable: its simple and elegant formulation [5] in 
terms of the multiple scattering theory [6] and the fact 
that, since it cleanly defines a homogeneous mean field 
alloy, it constitutes a natural starting point for pertur- 
bative studies of the fluctuations. The last feature has 
allowed the development of successful methods for the 
study of phase equilibria [7] and magnetic phenomena [8] 
in metallic alloys and led to surprisingly accurate calcu- 
lations of properties connected with typical Fermi liquid 
effects, such as spin [9] and concentration waves [10]. 

Since 1990, in spite of all the above successes, the CPA 
theory, or, better, the way in which it has been first im- 
plemented self-consistently within the density functional 
theory [11], has been criticised because it does not ac- 
count for the electrostatic energy that, in metallic alloys, 
arise from charge transfers [12]. Alternative models have 
been proposed able to cope with this lack [13] and to 
determine more accurately total energies and mixing en- 
thalpies. Unfortunately, the theories of this new class do 
not allow for an easy evaluation of spectral properties 
and are usually much more computationally demanding 
than the CPA. 

In the last ten years several attempts have been made 
to improve the CPA theory including electrostatic in- 
teractions. The starting point of these modified CPA 
theories is the consideration of the physical mechanism 



responsible of the screening of impurity charges in met- 
als, this is generally took into account introducing some 
screening length of the order of the nearest neighbours 
distance. Such theories, the screened impurity model by 
Abrikosov et al. [14] or the screened CPA method by 
Johnson and Pinski [15] or the more recent model by Za- 
harioudakis et al. [16], are able to achieve considerable 
improvements over the standard CPA, both for the total 
energies and the spectral properties [3,17]. 

In the last few years, the advent of order N electronic 
structure calculations [18,19], making feasible the study 
of large supercells containing hundreds of atoms, has led 
to a very remarkable result concerning the charge trans- 
fers in metallic alloys. Namely, Faulkner, Wang and 
Stocks [20,21] discovered that local charge excesses in 
random alloys are related by simple linear relations, the 
'qV laws, to the local electrostatic potentials. This new 
result, although not yet formally derived within the con- 
text of the density functional theory, has been confirmed 
by other order N calculations [3] and, indirectly, by pho- 
toemission experiments [22,23] and does not appear com- 
patible with previously proposed models [12,24]. It is 
interesting to observe that the existence of qV linear re- 
lations has been obtained also within a Thomas-Fermi 
model [25] and, henceforth, does not appear related with 
the specific form of the density functional used, but, 
rather, to the achievement of self-consistency in density 
functional theory calculations along with the solution of 
the Poisson equation used for the potential reconstruc- 
tion. 

The new results of Refs. [20,21] changed the scenario: 
the local charge excesses take all the possible values over 
some range, they can be described by a distribution that 
cannot be reproduced just by considering the nearest 
neighbours environment. The standard CPA theory or 
its modified versions [14-16], cannot describe such a dis- 
tribution, nor other kind of models [13] can. A possi- 
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ble way out is to renounce to the assumption of a single 
atomic potential for each alloying species, here referred as 
to the isomorphous assumption. Thus, for instance, the 
polymorphous CPA (PCPA) of Ujfalussy et al. [26,27], 
through the explicit inclusion of the electrostatic inter- 
actions between different crystal sites, leads to different 
site potentials. Though very successful from many points 
of view, however, the PCPA model, at similarity of order 
N calculations, is based on the use of supercells. As no- 
ticed by Faulkner [28], this does not forbid to obtain ac- 
curate ensemble averages for the relevant physical observ- 
ables, nevertheless the kind of averages that can be ob- 
tained from supercells is bound to some specified model 
of disorder, usually the random alloy model in which no 
correlation is assumed for the occupations of two differ- 
ent crystal sites. This, for instance, forbids the study of 
order-disorder phase transitions. 

In this paper we shall follow an alternative route. We 
shall analyse the response of the mean field CPA alloy 
to external local fields designed to simulate the effects of 
local electrostatic fields. Our main result will be that the 
qV relations can be obtained within the context of the 
isomorphous CPA theory, in satisfactory agreement with 
order N calculations. The theoretical scheme we shall 
develop allows for full charge relaxation at the impurity 
sites and clarifies that the linearity of the qV relations 
holds well beyond the linearity regime for other site prop- 
erties. Of course, in virtue of the fluctuation-dissipation 
theorem, the response to the external field that can be 
obtained from our theory must be equal to the response 
to the internal field due to the electrostatic interactions. 
Thus, our results could constitute a first step towards 
a CPA based theory that at the same time will be able 
to account accurately for the electrostatic interactions 
in metallic alloys and that maintain all the traditional 
strengths of the CPA. Another important feature of our 
work is that it does not rely on some specified supercell 
and, thus, in principle, can be used for studying ordering 
phenomena in metallic alloys. 

The remainder of this paper is organised as follows. In 
the next Section II we shall analyse the information about 
charge transfers that can be obtained from the standard 
isomorphous CPA theory in the light of the new results 
from order N calculations. In Section III we shall de- 
scribe a new version of the Coherent Potential Approx- 
imation that includes local fields (CPA+LF) as well as 
the applications of the model to fee CuPd and fee and 
bec CuZn alloys. In the same Section III we shall discuss 
the picture obtained for the charge relaxation phenom- 
ena, connected with the screening of the external field. 
In the final Section IV we shall make our summary and 
draw our conclusions. 



II. CHARGE TRANSFERS IN METALLIC 
ALLOYS. 

A. Charge transfers in the CPA theory. 

As noticed by Magri et al. [12], the Coherent Potential 
Approximation [2] does not include the energetic contri- 
butions arising from the charge transfers between differ- 
ent sites in metallic alloys [29]. This notwithstanding, 
sometimes, the information contained in the CPA the- 
ory about the same charge transfers, has been helpful 
for understanding the physical properties of certain sys- 
tems [30,31]. Below we shall elaborate about the reasons 
for this apparent paradox. 

The CPA theory (we use the multiple scattering the- 
ory formalism [5,32]) deals with a random binary alloy 
A CA B CB by solving for tc the so called CPA equation, 

c A G A {t A ,t c ) + c B G B {t B ,tc) = G c (t c ) (1) 

The three Green's functions in Eq. (1), G A (t A ,tc), 
GB(tB,tc) an d Gc(tc), refer to the three different 
problems sketched in Fig. 1. In fact, Gc(tc) is the 
Green's function for an infinite crystal whose sites are 
all occupied by effective scatterers characterised by the 
single-site scattering matrix tc- On the other hand, 
GA(B)(tA(B)i tc) is the Green's function for a single im- 
purity 'atom' described by the single-site scattering ma- 
trix t A ( B ) and embedded in an infinite crystal with all 
the other sites characterised by the single-site scattering 
matrix tc- While the homogeneous effective crystal, the 
'coherent' medium of the CPA theory, let us call it C, is 
electroneutral, the two impurity Green's functions lead 
to net charge excesses, q A and q B , in the sites occupied 
by the A and B impurities. On behalf of Eq.(l), these 
charge excesses satisfy the condition, 

c A q° A + c B q°B = (2) 

However, in the coherent medium C, there is no charge 
transfer from one site to the others, thus, Eq.(2) cannot 
be interpreted as an ordinary clectroneutrality condition. 
This notwithstanding, the excess charges, q A ( B y can le- 
gitimately be considered as the amount of charge that 
the impurity site A(B) would attract from the medium 
C [29], as in the two configurations on the LHS of Fig. 1. 
We could think that Eq. (2) describes an indirect charge 
transfer, from A to C and from C to B. Since C is a mean 
field approximation that accurately accounts for the elec- 
tronic properties of the alloy considered [3] , we could say 
that the CPA 'charge transfers' in Eq.(2) reflect the elec- 
tronegativity differences between A and B atoms in the 
random alloy A CA B CB and, hence, contain valuable physi- 
cal information. To put it into other words: in the physics 
of metallic alloys the mean field metal C plays the same 
role of reference system that, for molecules, is played by 
the H atom according to Pauling's celebrated concept of 
electronegativity [33]. A further reason of interest for 
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the CPA 'charge transfers', q A my ^ s * ne ^ ac ^ that these, 
though obtained from an isomorphous model, at similar- 
ity of the 'true' charge transfers, show up dependencies 
on the alloy concentration that cannot be reproduced by 
trivial models [30]. 

An additional difficulty in interpreting q A , B \ as the 
'true' charge transfers in a random alloy is that, on a 
more physical ground, the presence of a net charge on A 
should induce charges on the effective scatterers C and 
these, in turn, would affect the charge on A. Of course, 
as it is implied by its mean field nature, the CPA the- 
ory misses these important charge polarisation effects. It 
is easy to see that, the inclusion of these effects leads 
necessarily to a picture in which different sites occupied 
by atoms of the same kind are no longer equivalent be- 
cause they are affected in different ways by the charge 
rearrangements. In such a case, the isomorphous ap- 
proximation is no longer feasible. This circumstance led 
Ujfalussy et al. [26] to develop a 'polymorphous' coherent 
potential approximation (PCPA), in which different site 
potentials are allowed for. 

In this work we shall follow an approach alternative 
yet complementary to that of Ref. [26]. Borrowing the 
terminology by Ujfalussy et al., we shall include in the 
isomorphous CPA local perturbation fields and study the 
response of the charge density to these perturbations. 

B. Charge transfers in 'exact' multiple scattering 
theory calculations. 

In the last few years, Faulkner, Wang and Stocks 
[20,21] have extensively studied the distribution of 
charges that results from their Locally Self-consistent 
Multiple Scattering (LSMS) [18] calculations in binary 
metallic alloys. LSMS calculations are basically exact, 
except for the Local Density and muffin-tin approxima- 
tions used, and deal with large supercells containing 100 
to 1000 atoms, designed to simulate substitutional dis- 
order. The principle result of Ref. [20,21] is that, for a 
given supercell configuration, a simple linear law relates 
the net site charges, qi, to the Madelung potentials, Vi, 
at the same sites, say, 

a%qi + Vi = ki (3) 

For a binary alloy the quantities a* and ki take the values 
a a and Ua if the i-th site is occupied by a A atom or as 
and fee if it is occupied by B. The Madelung potentials 
Vi depend on the charges at all the crystal sites, through 
the relationship 

V i = 2j2M ij q j (4) 

i 

where the factor 2 comes from using atomic units. The 
Madelung matrix elements, My, are defined [34] as 



where the r;j are the translations from the i-th to the 
j-th site within the supercell and R are the superlattice 
translation vectors. 

Very remarkably, no numerically significant deviations 
from the linearity have been found, at least in the range 
of the variations of qi and Vi likely to occur in metallic al- 
loys. This result, as stressed by the same authors [20,21], 
is by no means obvious and should be understood as 
a consequence of the complex charge relaxation phe- 
nomena that occur, in the calculations, along with the 
achievement of self-consistency, or, in the physical world, 
through the mechanism mentioned in the previous sub- 
section. Actually, it has been shown that simple models 
assuming the site charges proportional to the number 
of unlike neighbours [12] lead to site charges distribu- 
tions [24] not compatible with Eq. (3). 

Although the above mentioned qV linear laws have 
been extracted from first principle calculations, their for- 
mal derivation within the density functional theory has 
not yet been obtained. Therefore, Eq. (3) should be con- 
sidered as an empirical law that holds, at least, within 
the basic approximation underlying the LSMS theory, 
i.e. the local density approximation for the density func- 
tional theory and the muffin-tin approximation for the 
crystal potential. On the experimental side, the validity 
of Eq. (3) is supported by photoemission spectroscopy 
measurements on random alloys [22,23]. 

An aspect of the above phenomenology, crucial for our 
present concerns, is the fact that, for random alloys, the 
four constants in Eq. (3) depend only on the site occupa- 
tions and the supercell configuration [35]. Moreover the 
same constants exhibit only minor variations for differ- 
ent configurations corresponding to the same mean alloy 
concentration. Remarkably, the dependence of the above 
constants on the molar fractions does not appear easy 
to describe within simple models. All the above circum- 
stances suggest us that the constants in Eq. (3), that 
appear to depend mainly on the site chemical occupa- 
tion, could probably be calculated within a isomorphous 
theory based on the CPA. Below we shall develop such 
an approach. 

III. RESPONSE TO LOCAL FIELDS OF THE 
'CPA ALLOY'. 

A. The local field CPA (CPA+LF) model 

As discussed in the previous section, LSMS calcula- 
tions suggest that the charge excess at the i-th site of a 
metallic alloy is determined, at least within the accuracy 
permitted by numerical errors, by the local Madelung 
field, Vi, and by the site occupation, say A or B, for a 
given alloy concentration. Borrowing the language of the 
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Ginzburg-Landau theory, we can think of the whole set 
of local charges, Qi, as the order parameter field. The 
phenomenology suggests a local view in which the local 
excess of charge, (ft, can be considered as the response to 
the local perturbation field, Vi. 

In this section we develop a simple model for the charge 
response to such local perturbations. At variance of what 
happens in real systems, for which the local field at the 
i-th site is determined by the charges at all the other sites 
(see, e.g. Eq. (4)), we shall treat the external local field, 
let us call it as a parameter that can be varied at will. 
We imagine to have an A impurity atom in a otherwise 
homogeneous crystal with all the other sites occupied by 
C scatterers. We suppose that the single site scattering 
matrix of the CPA medium, tc, and its Fermi energy, 
Ep, have been determined by the CPA theory for the 
binary alloy A CA B CB . The local external field, $, takes 
a constant value within the impurity site volume and is 
zero elsewhere. This situation is pictorially represented 
in Fig. 2. To simplify our discussion we shall solve the 
problem using the Atomic Sphere Approximation (ASA) . 
However, the following considerations hold for any cellu- 
lar method, and, with minor modifications, also for the 
muffin-tin approximation. 

We shall refer to the impurity A in the presence of 
the external field $ as to (A, $). When $ = 0, the site 
Green's function associated with it, G\(t\,tc), reduces 
to the usual CPA Green's function, GA{tA,tc)- When 
$ ^ 0, G\{t\,tc) can be readily obtained using the mul- 
tiple scattering theory impurity formula [32]: 

G* A {E,r,r>) = Y,[Zl{E,f)rl LL ,Zt,{E,r') 

L,L' 

Zt(E,r)J*(E,r')5 LL ,} (6) 

where 

T% = D*T C =[1 + T C ((t*)- 1 &)] ^ T C (7) 

In Eqs. (6) and (7), E is the energy, tc and tc are the 
CPA single site scattering matrix and scattering-path op- 
erator, as determined by an isomorphous CPA calcula- 
tion, i.e. $ = 0, for the alloy at hand. The single site 
scattering matrix corresponding to (A, $), £*, is to be 
determined from the site potential V®(r) + $, D\ is 
the CPA projector relative to the same site potential, 
Zf(E,r) and jf i {E,r) are two orthogonal solutions of 
the Schroedinger equation for the same potential, chosen 
in such a way that the first behaves regularly at r = 0. 
In our notation L = (I, m) labels the angular momentum 
quantum numbers and, for sake of simplicity, we omit the 
energy dependence of all the scattering matrices. A com- 
plete account of the notation can be found in Refs. [32] . 

The charge density corresponding to (A, $) is obtained 
integrating Eq. (6) over the energy to the Fermi level, 

p\{r) = -Urn | J** dE G\{E, r, f ' = f) j (8) 



The corresponding site potential, V®(r), can be recon- 
structed by solving the appropriate Poisson equation 
and adding the exchange-correlation contribution [36,37]. 
Unless $ = 0, it will be different from the site poten- 
tial obtained from the zero field CPA theory, V A (r) = 
Vjf =0 (r) , due to charge relaxations expected to screen in 
part the external field. In a numerical implementation 
of the theory, Eqs. (6-8) and the potential reconstruc- 
tion need to be iterated starting from a convenient initial 
guess, until convergence is achieved for Vjf (r) or, equiv- 
alently, for p\{r). Hereafter we shall refer to the above 
model as to the Local Field CPA (CPA+LF). 

Once convergence is obtained for the charge density, 
the net charge on the site A can be obtained by integrat- 
ing over the atomic sphere volume and subtracting the 
nuclear charge, Z A , 

q A (*) = J drp\(f)-Z A (9) 

It is important to realise that, while the above self- 
consistent determination of V% (r) or p\ (r) allows for full 
charge relaxation at the impurity site, the CPA+LF does 
not modify the properties of the CPA medium C: these 
remain specified by the quantities tc and Ep determined 
at zero external field. In other words, the charge relax- 
ation that is likely to occur also in the neighbouring sites 
is not accounted for. 

Although in this paper we deal only with random al- 
loys, we like to mention that our formalism can be ap- 
plied also to the study of the screening phenomena in 
pure systems. In this case, of course, the CPA scattering 
matrices, tc and tc, in Eq. (7) should be replaced by 
their counterparts for a pure system. 

B. CPA+LF results for CuZn and CuPd alloys: the 
charge relaxation 

We have implemented the CPA+LF theory within our 
well tested KKR-CPA code [38]. If t c and t c from a 
previous standard KKR-CPA calculation are stored on a 
convenient energy mesh, the extra computational efforts 
required by the CPA+LF calculation are negligible. 

In this paper we discuss results for fee CuPd and for 
bec and fee CuZn random alloys at several concentra- 
tions. In all the cases we have used the Local Density 
approximation (LDA) for the exchange-correlation po- 
tential [39], the ASA approximation for the site poten- 
tials and the angular momentum expansions have been 
truncated at I max = 3. We have a fully relativistic 
treatment for core electrons and use a scalar relativistic 
approximation for valence electrons. For all the alloy sys- 
tems considered in this paper, the lattice parameters have 
been kept fixed on varying the concentration. In partic- 
ular, we set a = 5.5 a.u. and a = 6.9 a.u. for bec and 
fee CuZn, and a = 7.1 a.u. for fee CuPd. The reasons 
for this choice are, in first place, avoiding the considera- 
tion of the lattice parameter variations, not relevant for 
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the purpose of this paper, and, in second place, allowing 
for an easier comparison with the results of Refs. [20,21]. 
For future reference, we mention that with this choice, 
the atomic volumes in fee and bec CuZn alloys differ only 
about 1.3 per cent. 

As we said, the CPA+LF model allows for the deter- 
mination of the response to an external potential field by 
the electrons inside the atomic sphere A. More specifi- 
cally, the difference 

AVt(r) = V*{r) + $ - V*=°(r) (10) 

can be interpreted as the sum of the external field, $, 
and the internal screening field inside the atomic sphere. 
Some typical trends for this quantity are shown in Fig. 3. 
There we report AV c f(r), (a — Cu,Pd), for an fee 
Cuo.50Pdo.50 random alloy, that we have selected as a 
typical case. At the Wigner-Seitz radius, rws ~ 2.77 
a.u., the internal field is able to screen about one half of 
the external field, both for Cu and Pd impurities, while 
the screening is almost complete for r < 1 a.u.. Appar- 
ently, the effect of the screening is far from being just 
a constant shift of the local chemical potential: if that 
was the case, in Fig. 3 we would have just equally spaced 
horizontal lines. What we observe is much more com- 
plicated. For instance, in the case of Pd impurities we 
see that, for the largest values of $ considered, AVp d (r) 
takes negative values at small r's, thus, in the case of Pd, 
the external field appears overscreened in the same range 
of r. 

The complex nature of the screening phenomena is fur- 
ther confirmed by a look at the electronic densities. In 
Fig. 4 we plot the excess charge density induced by the 
external field 

A/>S(r)=p$(r)-pr (0 (11) 

both for Cu and Pd sites, again for random fee 
Cuo.50Pdo.50. The largest effects come from the large 
r region, where the electron density decreases on increas- 
ing $ (everywhere in this paper the expressions "elec- 
tronic density" or " charge density" are used indifferently 
with the meaning of "electron number density", i.e. the 
charge factor, — e, is not included). In the innermost part 
of the atomic spheres, the variations of the charge den- 
sity sometimes may have opposite sign with respect to 
that observed close to the cell boundary. 
We have considered also the quantity, 

«M - pt{ t*X) r) ™ ^ l09ptir) (12) 

that, in the limit $ — > reduces to the logarithmic 
derivative of p*(r) and that, on the basis of a formal 
scattering theory analysis [40] is expected to have a weak 
dependence on $. As we can see from Fig. 5, the residual 
dependence on $ is about a few per cent in a relatively 
small r interval not far from r w $ an d less then 1 per 



cent in most of the atomic sphere. Although the infor- 
mation conveyed by Fig. 5 can be valuable for the purpose 
of improving the initial guesses for the charge densities, 
however the dependence of 6* (r) on r appears still quite 
complicated: 6*(r) is very small at small r's and takes 
its largest (negative) value at about 2.3 Bohr radii for Cu 
and 2.4 Bohr radii for Pd. In this region of the atomic 
sphere the charge density appears most sensitive to the 
external field. 



C. CPA+LF results for CuZn and CuPd alloys: the 
site charges 

As discussed in the previous subsection, the charge re- 
laxation phenomena occurring when the local field $ is 
included in the isomorphous CPA theory are quite com- 
plicated. It is then surprising to see that the correspond- 
ing site charges, calculated from Eq. 9, exhibit much sim- 
pler trends. 

In Fig. 6 we report q a (a =Cu, Pd) vs. $ for a 
Cuo.50Pdo.50 fee random alloy. As it is evident, the data 
can be fitted very well by two straight lines, one for each 
atomic species. The resulting correlations differ from one 
by less than one part over a million. Interestingly, the 
slopes of the two lines arc different by a relatively small 
but statistically relevant amount, about 5 per cent. We 
notice that in Fig. 6 we have considered also <j> values 
considerably larger that those observed in LSMS calcu- 
lation or likely to occur in real systems, thus our data 
support the view that the linearity observed has little to 
do with the size of $. 

We have fitted the q a vs. $ curves at each molar frac- 
tion for fee Cu c Pdi_ c , fee Cu c Zni_ c and bec Cu c Pdi_ c 
random alloys, at a number of alloy concentrations, using 
the linear relationships 

g Q ($) =q° a -R a $ (13) 

However, at $ = 0, our CPA+LF model satisfies the CPA 
'electronegativity' condition, Eq. (2). Henceforth, q A and 
q B are not independent quantities and we have chosen as 
the parameters of our fit only the three quantities Ra, 
Rb and 

A = q° A /c B = - q B /c A = q A -q B (14) 

The results of these fits are reported in Table I. The de- 
viations from linearity are always comparable with the 
numerical errors and the fits cannot be improved signifi- 
cantly using different functional forms. The trends found 
for the fitting parameters vs. the alloy molar fractions are 
shown in Fig. 7. The dependence on the concentration 
is appreciable for all the fitting parameters, as expected 
on the basis of the arguments developed in Section II. 
Remarkably, the dependences on the alloy system and 
on the concentration appear at least as much important 
as that on the atomic species. Thus, for instance, for a 
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given alloy system and concentration, there are relatively 
small differences between the values of R corresponding 
to sites occupied by different atoms. On the other hand, 
we find much larger variations for Rc u throughout the 
alloy systems considered. It is interesting to observe that 
the trends for the slopes, Rc u and Rz n , and for A are 
very similar in both fee and bec CuZn alloys. We notice 
also that A, a measure of the electronegativity difference 
between the alloying species, exhibits, at least for CuPd 
alloys, non negligible variations vs. the concentration. 
In the model of Ref. [12], the same quantity is assumed 
independent on the concentration. As we see from Table 
I, the values for A from our theory are systematically 
25-35 per cent smaller than those from LSMS calcula- 
tions. This is, actually, a feature of the standard CPA 
theory, because the local fields do not enter in the deter- 
mination of A and it has already been discussed in the 
literature [31]. This notwithstanding, the CPA is able to 
catch the qualitative trends of A vs. the concentration 
for all the systems considered. 

The LCPA+LF model gives for q vs. $ the same linear 
functional form as that obtained for q vs. V from LSMS 
calculations. However, the differences between the two 
different sets of calculations forbid, at this stage, a di- 
rect comparison of the fit coefficients. In fact, as we have 
already stressed, our CPA+LF model does not account 
for charge relaxation outside the impurity site volume. 
By its construction, the CPA medium C is able to screen 
the impurity charge at $ = 0, i.e. We can think that 
this amount of charge is screened by the infinite volume of 
C. The introduction of the local field at the impurity site 
causes a local excess of charge, q a (&) — g° , with respect to 
the standard CPA. In order to have global electroneutral- 
ity in the CPA+LF theory, we have to introduce, some- 
where outside the impurity site, an opposite amount of 
charge, — q a ($>). This can be accomplished follow- 
ing the ideas contained in the screened impurity model 
(SIM-CPA) by Abrikosov et al. [14] and supposing that 
the excess (with respect to the standard CPA) charge at 
the impurity site, q a ($) — 3°' is completely screened at 
some distance, p, of the order the nearest neighbours dis- 
tance, r\. Accordingly, each of the n nearest neighbours 
of the impurity cell shall have, in the mean, a net charge 
excess (g° — q a (&))/n. This, in turn, will cause an extra 
field $! = n(2/p)(q° a - « a ($))/n = 2(g° - g„($))/p on 
the impurity site. The total field at the impurity site will 
be then the sum of the external field $ and of the above 
extra term, in formulae, 

V r a =* + 2(g£-«a(*))/p (15) 

Then, by solving for $ the last equation and substituting 
in Eq. (13), we find 

qa{<P) = q ° a ~ l + 2R a /p Va = q «- RaVa (16) 

The coefficients R a can be compared directly with the 
slopes of the qV relations from LSMS calculations. How- 
ever, the comparison, reported in Table I, requires a 



caveat: we have assumed p = n, i.e. a complete screen- 
ing at the distance of the nearest neighbours. Actually, 
the screening lengths in metals are of the order of this 
distance [41], but our estimate is too crude to expect for 
a very good quantitative agreement with LSMS calcula- 
tions in which the charge relaxation is allowed at all the 
length scales. However, the agreement found is quite sat- 
isfactory, within 10 per cent, for CuPd alloys, while larger 
discrepancies are found for CuZn. Again, the trends for 
R a vs. the concentration are qualitatively reproduced. 

IV. CONCLUSIONS 

Using a convenient generalisation of the CPA theory 
to include a local external potential, we have been able 
to reproduce the peculiar linear relationship between the 
local charge excesses and the local Madelung fields, ob- 
tained previously from the analysis of LSMS data [20,21]. 
Our CPA+LF theory has a very simple structure and, 
as opposite to LSMS, requires modest computational ef- 
forts. The comparison of the linear laws coefficients with 
those from LSMS calculations is quantitatively satisfac- 
tory, though there is place for further improvement. In 
particular, it has been necessary to include some screen- 
ing length, that we have kept fixed to the nearest neigh- 
bours distance, and that constitute the only non ab ini- 
tio input of our theory [3]. Work is in progress about 
this point, and we hope to present in the next future a 
fully ab initio version of the theory. We think that the 
CPA+LF model catches the basic physical mechanisms 
that determine the distribution of the charge excesses in 
metallic alloys and we are currently working on a new 
version of the theory in which the distribution of local 
charge excesses and that of the local fields are obtained 
self-consistently. We are confident that this will allow 
for an accurate calculation of the electrostatic contribu- 
tion to the total alloy energy and believe that these new 
techniques can prove useful for the study of the phase 
equilibria in metallic alloys. 

The main result of the present paper, having obtained 
the correct functional form for the qV relations within a 
CPA-based theory, does not solve the problem of their 
derivation within the density functional theory Never- 
theless, our work offers a much more comfortable math- 
ematical ground for the search of such a solution and 
clarifies two points that we like to list below. 

i) A conceptual advantage of our CPA+LF theory is 
that it is build upon a cleanly defined reference medium, 
the CPA alloy, that is kept fixed. Having fixed the refer- 
ence medium, the source of the linear laws can be traced 
only in the CPA projectors and in the site wavefunctions 
(see Eqs. (6) and (7)), these, in turn, are determined by 
the nuclear charge on the impurity site and the coupling 
potential entering in the corresponding Schroedinger- 
Kohn-Sham equation: in other words, by the chemical 
species and the Madelung potential. This ensures that, 
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for a given alloy, i.e. for some specified alloy Green's 
function, any site property depends only on the chem- 
ical species and the Madelung potentials. This applies 
also, to some extent, to the more realistic LSMS model 
in which, having fixed the alloy configuration, fixes the 
system's Green's function. Of course, in the LSMS model 
the Green's function is determined self-consistently, while 
in the CPA+LF model we use a reasonable approxima- 
tion of it: this difference, however is not essential since in 
both cases there is a unique Green's function. A more im- 
portant difference between CPA+LF and LSMS calcula- 
tion is the fact that, in the first theory, the way the impu- 
rity site is embedded in the system is determined simply 
by the CPA projectors, while in a more complete treat- 
ment of the multiple scattering, as in the case of LSMS 
calculations, the relationships between the site-diagonal 
part of the Green's function at the impurity site and the 
full system's Green's function are more complicate. Also 
in the last case one could try to define something analo- 
gous to the CPA projectors, such generalised projectors, 
however will depend not only on the systems Green's 
function and on the impurity site properties but also on 
the environment of the neighbouring sites. The fact that 
the qV laws have been obtained first from the analysis of 
LSMS data suggests that corrections for the neighbour- 
ing sites are really important only for the electrostatic 
coupling term. This view is supported also by the very 
remarkable accuracy achieved by the polymorphous CPA 
theory [26]. 

ii) Provided that the effect of the nearest neighbours 
is negligible, on the basis of the above point i), any site 
property is a unique function of the local field and chem- 
ical occupation. This, together with the linear response 
theory, of course, implies linear qV relations for small 
fields. However, the range of fields investigated in this 
paper has been sufficiently wide to encounter non lin- 
ear behaviours for many quantities, as, for instance, the 
value of the charge density at some specified r. Although 
the reasons for this remain elusive, we can conclude that 
the trends for q remains linear beyond the linear response 
regime for p(r). 

After the submission of the present work, two pa- 
pers [42,43] appeared that treat some of the questions 
addressed here from a different point of view. In Ref. [42] 
it is demonstrated that linear qV laws can be obtained 
within the SIM-CPA model of Abrikosov et al. [14] on 
varying the effective screening length contained in the 
same model. Since the SIM-CPA assumes an effective 
field proportional to the charge excess, and the pro- 
portionality constant is the inverse of the above effec- 
tive length, it is clear that the procedure described in 
Ref. [42] has effects similar to varying a local external 
field. Though we believe that the results of Refs. [42,43] 
offer an independent confirmation of ours, there are im- 
portant differences between the two approaches. To men- 
tion but one: just because of its mathematical construc- 
tion, the SIM-CPA model forces the qV laws to have 
the same slopes for both the alloying species, while the 



LSMS calculations of Ref. [20,21], in agreement with the 
CPA+LF calculations presented in this paper, find dif- 
ferent slopes for different alloying species. 
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Alloys 


c 


A 


Rcu 


Rx 


RMS xlO 4 


Rcu 


Rx 


A 


Rcu 


Rx 


fee Cu c Pdi_ c 


0.10 


0.183 


1.093 


1.156 


1.8 


0.762 


0.792 


0.238 


0.833 


0.843 




0.25 


0.175 


1.124 


1.187 


2.1 


0.776 


0.806 


0.229 


0.838 


0.851 




0.50 


0.160 


1.184 


1.244 


1.9 


0.805 


0.832 


0.219 


0.843 


0.851 




0.75 


0.150 


1.243 


1.288 


2.4 


0.831 


0.851 


0.212 


0.838 


0.853 




0.90 


0.148 


1.267 


1.307 


4.4 


0.842 


0.860 


0.211 


0.836 


0.853 


bee Cu c Zni-c 


0.10 


0.109 


1.206 


1.232 


10 


0.800 


0.812 


0.155 


0.536 


0.581 




0.25 


0.114 


1.237 


1.255 


10 


0.814 


0.822 


0.159 


0.526 


0.554 




0.50 


0.116 


1.237 


1.251 


6.9 


0.814 


0.820 


0.156 


0.545 


0.549 




0.75 


0.116 


1.247 


1.255 


5.0 


0.819 


0.822 


0.155 


0.567 


0.564 




0.90 


0.116 


1.248 


1.254 


3.2 


0.819 


0.822 


0.158 


0.582 


0.577 


fee Cu c Zni_ c 


0.10 


0.106 


1.202 


1.223 


8.2 


0.805 


0.815 


0.145 


0.575 


0.628 




0.25 


0.111 


1.220 


1.237 


8.1 


0.813 


0.821 


0.150 


0.580 


0.618 




0.50 


0.116 


1.222 


1.241 


5.5 


0.814 


0.822 


0.151 


0.600 


0.622 




0.75 


0.117 


1.247 


1.256 


5.2 


0.825 


0.829 


0.150 


0.615 


0.632 




0.90 


0.118 


1.249 


1.256 


3.3 


0.826 


0.829 


0.152 


0.616 


0.630 



TABLE I. Fit parameters for the q vs. $ relationships from CPA+LF calculations in fee Cu c Pdi_ c , bee Cu c Zni- c and fee 
Cu c Zni_ c random alloys. The 'electronegativity difference', A, and the response coefficients, Rcu and Rx, {X — Pd or Zn, as 
convenient) are defined in Eqs. (13) and (15), RMS is the root mean square deviation. The 'renormalised' response coefficients, 
Rcu and Rx, are defined in Eq. (17). On the right "we report A, Rcu and Rx from the LSMS calculations of Refs. [19-20]. 
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FIG. 1. Schematic illustration of the CPA theory. Dark 
sites are occupied by the CPA coherent scatterer described 
by the single site scattering matrix tc- The central impurity 
sites, labelled by A and B, are characterised by the single site 
matrices tA and is- 



FIG. 2. Schematic illustration of the CPA+LF method. As 
in Fig. 1, dark sites are occupied by the CPA coherent scat- 
terer described by tc- In the central site, occupied by A, acts 
also a constant field $. 
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FIG. 3. Calculated total field AV*(r),a =Cu, Pd (see Eq. 
10) in Cuo.50Pdo.50 fee random alloys. The labels indicate the 
values of the external field, <E>. At the Wigner-Seitz radius, 
rws ~ 2.77 a.u., the total field results to be about one half 
of the external field, while the electronic screening is almost 
complete at r — 0. 
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FIG. 4. Calculated excess charge density, Airr 2 Ap®(r) 
(a =Cu, Pd) (see Eq. 11) in Cuo.50Pdo.50 fee random alloys. 
The labels indicate the values of the external field, $. 
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FIG. 5. The 'logarithmic derivative' with respect to the 
external field(Eq.(12)), b*(r) (a=Cu,Pd) in Cuo.50Pdo.50 fee 
random alloys. The continuous and the dotted lines refer, 
respectively, to $ = —0.3 and $ = 0.3, i.e. the lowest and the 
highest <E> values considered in Fig.3. 
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FIG. 6. Site charge excesses q a (a = Cu,Pd) vs. the exter- 
nal field, from CPA+LF calculations for Cuo.50Pdo.50 fee 
random alloys. Circles and diamonds, respectively, indicate 
Cu and Pd impurities. 
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FIG. 7. Fit coefficients of the linear law q vs. <!> from 
CPA+LF calculations for fee CuPd and fee and bcc CuZn 
alloys plotted vs. the Cu content. Upper frame: response co- 
efficients R a , (a refers to the alfoying species); iower frame: 
'efectronegativity difference', A. The various alioy systems 
are indicated by iabefs. 
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